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Abstract 

We present a range of physics results for central exclusive production processes at 
the LHC, using the new SuperChic 2 Monte Carlo event generator. This includes 
significant theoretical improvements and updates, most importantly a fully differential 
treatment of the soft survival factor, as well as a greater number of generated processes. 
We provide an overview of the latest theoretical framework, and consider in detail a 
selection of final states, namely exclusive 2 and 3 jets, photoproduced vector mesons, 
two~photon initiated muon and W boson pairs and heavy Xc,b quarkonia. 


1 Introduction 

Central Exclusive Production (GEP) is the reaction 

pp{p) p + X + p{p) , 

where ‘+’ signs are used to denote the presence of large rapidity gaps, separating the system X 
from the intact outgoing protons (anti-protons). Over the last decade there has been a steady 
rise of theoretical and experimental interest in studies of this process in high-energy hadronic 
collisions, see mu for reviews. Theoretically, the study of GEP requires the development of 
a framework which is quite different from that used to describe the inclusive processes more 
commonly considered at hadron colliders. Moreover, the dynamics of the GEP process leads 
to unique predictions and effects which are not seen in the inclusive mode. Experimentally, 
GEP represents a very clean signal, with just the object X and no other hadronic activity 
seen in the central detector (in the absence of pile up). 


The CEP process requires the t-channel exchange of a color-singlet object, so that the 
outgoing protons can remain intact. One possibility to achieve this is the two-photon fusion 
process 77 —)■ X, where the radiated photons couple to the electromagnetic charge of the 
whole protons. Alternatively, the process may be mediated purely by the strong interaction: 
provided the object X mass is large enough, this can be considered in the framework of 
pQCD, via the so-called Durham model is]. Finally it is possible for ‘photoproduction’ 
reactions to occur, where one proton interacts electromagnetically and one interacts strongly. 

In any detailed phenomenological study of such processes, it is important to have a Monte 
Carlo (MC) implementation, so that theoretical predictions can be compared more directly 
with experimental measurements. For this reason the authors have previously produced 
the publicly available SuperChic MC [5l[6], for the CEP of lighter Standard Model (SM) 
objects within the Durham model. While first considering Xc,b and Pcb quarkonia, this has 
subsequently been extended to include 77 and light meson pair {7ni,ri{')r]{')...) production, 
as well as the photoproduction of C-odd vector mesons. Such processes have been measured 
and compared to the MC predictions at the Tevatron and LHC, see for example |7HS], with 
results that are generally in good agreement. Other related available generators are the 
ExHuME [10] and FPMC [H] MCs. 

However, there exist a wider range of processes that are not included in earlier versions 
of SuperChic, but which have much phenomenological relevance, in particular in the light of 
the measurement possibilities for exclusive processes during Run-H of the LHC [T2|. Here, 
exclusive events may be measured with both protons tagged using the approved and installed 
AFP [13] and CT-PPS [TT] forward proton spectrometers, associated with the ATLAS and 
CMS central detectors, respectively, see also [15], as well as using rapidity gap vetoes to select 
a dominantly exclusive event sample. This latter possibility is in particular relevant at LHCb, 
for which the relatively low instantaneous luminosity and wide rapidity coverage allowed by 
the newly installed HERSCHEL forward detectors [H] are highly favourable, while similar 
scintillation counters are also installed at ALICE im. 

Particularly relevant is the case of exclusive jet production, which has been observed by 
both CDF [18] and DO [19] at the Tevatron, and for which already a sample of ‘exclusive-like’ 
2 and 3-jet events has been collected in a combined CMS+TOTEM run at 8 TeV [2U1I^ : 
exclusive 3-jet production in particular has not been implemented in any public MC. A 
further topical process is the exclusive production of quarkonia pairs, measured by LHCb 
in [ 22 ] and considered theoretically in [23], but for which a MC implementation has previously 
not been made publicly available. 

In addition to including a wider range of processes such as these, there are a number of 
theoretical updates and improvements which it is important to consider. Most signihcantly, in 
all exclusive processes it is necessary to account for the probability that no additional particles 
are produced by soft proton-proton interactions, independent of the hard process: the so- 
called ‘survival factor’. In previous versions of SuperChic as well as in other generators mm 
this is simply treated as a constant probability which suppresses the overall cross section. 
However, it is well known that the survival factor can depend sensitively on the hnal-state 
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particle momenta, and so such an averaging will omit the influence this can have on the 
predicted distributions as well as only providing an approximate estimate of the (process- 
dependent) overall suppression. 

With these considerations in mind, we present in this paper results of the new SuperChic 
2 MC generator. This contains a range of theoretical improvements compared to the previ¬ 
ous version, most significantly including a fully differential treatment of the survival factor, 
maintaining the explicit dependence of this on the particle momenta in all cases. As well 
as the processes generated in the original MC, exclusive 2 and 3 jet, quarkonia {J/i’ and 
'ip{2S)) pair, SM Higgs boson production and the photoproduction of p and 0 mesons are 
now implemented. In addition, the two-photon production of 77 , W^W~ and lepton pairs 
are included; this is the first MC implementation of such photon-induces processes which 
includes a complete treatment of soft survival effects. The case of a e’''e“ initial state is also 
implemented for these processes. 

The outline of this paper is as follows. In Section |2] we summarise the theoretical ingredi¬ 
ents of the Durham model of QCD-mediated CEP, and describe how the soft survival factor 
can be included differentially in theoretical predictions and in a MC. In Section [3] we describe 
the theory of photon-induced processes, again providing details of how a full treatment of 
the survival factor can be achieved. In Sections 14.11 to 14.41 we present results of this MC for 
a range of processes: exclusive 2 and 3 jet production in Section 14.1( exclusive vector meson 
photoproduction in Section 14.2| two-photon induced W^W~ and lepton pair production in 
Section 14.3( heavy Xc,b quarkonia production in Section 14.41 In Section 14.51 a summary is 
presented of all processes that are generated, including some motivation for specific measure¬ 
ments that may be performed at the LHC. In Section |5] we briefly describe the SuperChic 2 
MC and its public availability. Finally, in Section [ 6 ] we present a summary and outlook. 


2 QCD processes 

2.1 Basic formalism 


CEP processes that proceed purely by the strong interaction can be described by the ‘Durham’ 
model, a pQCD-based approach that may be applied when the object mass Mx is sufficiently 
high, see PEllM] for reviews. The formalism used to calculate the perturbative CEP cross 
section is explained in detail elsewhere [Tll^ [251 - [2^ and we will only present a very brief 
summary here. The perturbative CEP amplitude, corresponding to the diagram shown in 
Fig. [H can be written as 


T = TT^ 




Qi(Q± - Pu)^(Q± + P2j 


fg{xi,x[,Ql,pl;ti)fg{x2,X2,Ql,pl;t2) , ( 1 ) 


where Q± is the transverse momentum in the gluon loop, with the scale = Q\ in the 
forward proton limit (see e.g. | 6 ] for a prescription away from this limit), and Ai is the 
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Figure 1: The perturbative mechanism for the QCD-induced exclusive process pp ^ p + 
X + p, with the eikonal and enhanced survival factors shown symbolically. 


colour-averaged, normalised sub-amplitude for the gg ^ X process 


M = 


2 1 
M|. - 1 




( 2 ) 


Here a and b are colour indices, Mx is the central object mass, is the gg ^ X vertex, 
are the transverse momenta of the incoming gluons, and ti is the squared momentum 
transfer to the outgoing protons. The /^’s in ([T]) are the skewed unintegrated gluon densities 
of the proton. These correspond to the distribution of gluons in transverse momentum Q±, 
which are evolved in energy up to the hard scale fip, such that they are accompanied by 
no additional radiation, as is essential for exclusive production. While the gluon momentum 
fractions Xi are set by the mass and rapidity of the hnal state, the fractions x/ carried by 
the screening gluon must in general be integrated over at the amplitude level. However, for 
the dominant imaginary part of the amplitude we have x' <C x, and it can be shown that the 
/g’s may be simply written as 


fg{x, X , Q_|_, Pp) 


d 


d\n{Q] 


Hr. 




(3) 


where Hg is the generalised gluon PDF [30], which for CEP kinematics can be related to 
the conventional PDFs [29l[3l]. The Tg in (I3l) is a Sudakov factor, which corresponds to the 
probability of no extra parton emission from each fusing gluon. 

We can decompose ([2]) in terms of on-shell helicity amplitudes, neglecting small off-shell 


4 














corrections of order ^ Omitting colour indices for simplicity, this gives 


' -|(qu'q2r)(J’++ + J’-) (^f = o+) 

-tl(qu xq2j|(T++-T__) (Jf = O') 

+|((?u^2^ - CO + *(CC + CO)^-+ (C = +2+) 
c +|((CC - CO - *(CC + CO)^+- (C = -2+) 


where the C indicate the parity and spin projection on the gg axis, and Tx^x 2 nre the 
corresponding g{Xi)g{X 2 ) X helicity amplitndes, see [SIE] for more details. In the forward 
proton limit (i.e. with ontgoing proton p± = 0) the only non-vanishing term after the Q± 
integration ([1]) is the hrst one: this is the origin of the selection rnle [mi32ll33] which operates 
in this exclnsive process, and strongly favonrs C = 0^ qnantnm nnmbers for the centrally 
produced state. More generally, away from the exact forward limit the non-J^ = O’*" terms 
in (j3]) do not give completely vanishing contributions to the integral and we hnd that 

ICI^.|=2)P {pl? 

|g4(J, = 0)|2 ~ (Qi)2 ’ ^ ^ 

which is typically of order ~ 1/50 — 1/100, depending on such factors as the central object 
mass, c.m.s. energy xfs and choice of PDF set |3l|6]. The on-shell decomposition (jl]) is used 
throughout, unless otherwise stated. 


2.2 Soft survival effects 

The expression ([I]) corresponds to the amplitude for the exclusive production of an object 
X in a short-distance interaction, that is, with no further perturbative emission. However 
secondary particles may also be produced by additional soft proton-proton interactions, 
independent of the hard process. Such underlying event activity will spoil the exclusivity of 
the event, and the probability that no additional particles are produced by accompanying 
soft proton-proton interactions is encoded in the so-called ‘survival factor’, see e.g. [3^1135] 
for some more recent theoretical work, and [3] for further discussion and references. 

The survival factor is not a simple multiplicative constant [6], but rather depends quite 
sensitively on the outgoing proton transverse momenta. Physically, this is to be expected, 
as the survival factor will depend on the impact parameter of the colliding protons; loosely 
speaking, as the protons become more separated in impact parameter, we should expect there 
to be less additional particle production, and so for the survival factor to be closer to unity 
(consequently, as we will see below, the average survival factor is much larger in the case of 
photon-mediated processes, where larger impact parameters are favoured, when compared to 
purely QCD processes). As the transverse momenta of the scattered protons are nothing 
other than the Fourier conjugates of the proton impact parameters, hu, we therefore expect 
the survival factor to depend on these. 


5 







For this reason, survival effects are included fully differentially in the final-state momenta 
in SuperChic 2. To describe in more detail how this is achieved, we can consider a simplified 
‘one-channel’ model, which ignores any internal structure of the proton; see [36l|37] for 
discussion of how this can be generalised to the more realistic ‘mutli-channel’ case. The 
average suppression factor is written as 

/o2 \ /d^bitd^bat |T(s,bit,b2i)Pexp(-f](s,6t)) 

^ /d2biid2b2dT(s,bii,b2t)|2 ’ 

where bj^ is the impact parameter vector of proton i, so that b^ = bi^ + h 2 t corresponds 
to the transverse separation between the colliding protons, with bt = |bt|. T(s,bit,b 2 t) 
is the CEP amplitude ([1]) in impact parameter space, and f2(s, bt) is the proton opacity, 
which can be extracted from such hadronic observables as the elastic and total cross sections 
as well as, combined with some additional physical assumption about the composition of 
the proton, the single and double diffractive cross sections. From (|6]), we can see that 
physically exp(—ff(s, 6^)) represents the probability that no inelastic scattering occurs at 
impact parameter bt. 

In the expression above, T(s, bii, h 2 t) is just the Fourier conjugate of the CEP amplitude 
dD. i.e. we have 

r(s,pi,,P2j = J d2bitd2b2te*P^x'bi*e-P2,-b2.y(^^bit,b2i) . (7) 


In transverse momentum space, the CEP amplitude including rescattering effects, T’’®®, is 
calculated by integrating over the transverse momentum carried round the Pomeron loop 
(represented by the grey oval labeled in Fig. [1]). The amplitude including rescattering 
corrections is given by 


lS,PU,P2xj = 


d^ki 

Stt^ 


Tei(s,ki) T(s, p'. ,p '2 


( 8 ) 


where p\^ = (pi^^ — k_i_) and p' 2 j^ = (p 2 j^ + k_i_), while T®\s, k^) is the elastic pp scattering 
amplitude in transverse momentum space, which is related to the proton opacity via 


r.i(s,() = 2s / d"b,e“‘ ''T.i(s,(),) =2« / d=b,e“‘ '=‘ 


(9) 


where t = —k^. We must add (|8]) to the ‘bare’ amplitude excluding rescattering effects to 
give the full amplitude, which we can square to give the CEP cross section including eikonal 
survival effects 


dcT 

d2pi^d2p2^ 


oc |T(s,pi^ 


P2j 


(S,Pi^,P 2 


( 10 ) 


where here (and above) we have omitted the dependence of the cross section on all other 
kinematic variables for simplicity. The overall normalisation of the cross section is achieved 
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exactly as in the unscreened case. It is this expression, suitably generalised to the multi¬ 
channel case, which is used in the MC. We note that following the discussion above, the 
expected soft suppression can be written in transverse momentum space as 

/<;2 \ _ |r(s.Pu.P2i)+r'“(s,Pu,P2i)P , , 

' /d2pi,d2p2jT(i,,pi,,P2jP ■ ^ ' 

It can readily be shown that ([6]) and (fTTD are equivalent. As expected, the soft suppression 
factor depends on the proton transverse momenta, and so may have an important effect on 
the distributions of the outgoing proton momenta, via (fTOj) . see e.g. [381139] . 

Besides the effect of eikonal screening S'eik, there is some suppression caused by the rescat¬ 
terings of the protons with the intermediate partons [331101111] • This effect is described by 
the so-called enhanced Reggeon diagrams and usually denoted as see Fig. [T1 The pre¬ 
cise size of this effect is uncertain, but due to the relatively large transverse momentum (and 
so smaller absorptive cross section of the intermediate partons, it is only expected to 
reduce the corresponding CEP cross section by a factor of at most a ‘few’, that is a much 
weaker suppression than in the case of the eikonal survival factor. Due to this uncertainty, 
in the current version of the MC these effects are omitted entirely, however by observing 
any departure from the MC predictions, for example in the invariant mass Mx distributions, 
such enhanced survival effects may still be investigated. 

Finally, it is worth pointing out that such soft survival effects do not only manifest 
themselves in CEP reactions. For example, the eikonal model of absorption discussed above 
may be used to predict the behaviour of the leading neutron spectra in diffractive dijet 
photoproduction, measured using a leading proton detector at HERA, see e.g. [^2]j33] . In 
addition, this approach may be used to explain [35] the breaking of factorisation for diffractive 
dijet production at the Tevatron [16], when this is given in terms of diffractive structure 
functions measured at HERA. As we only consider CEP reactions in this paper, we will not 
discuss these possibilities further here. 


3 Photon mediated processes 

3.1 Basic formalism 


Exclusive photon-exchange processes in pp collisions are described in terms of the equivalent 
photon approximation mi. The quasi-real photons are emitted by the incoming proton 
z = 1, 2 with a flux given by 


diV^(6) 


a d^^ 

TT + i'fml 









( 12 ) 


where and qi^ are the longitudinal momentum fraction and transverse momentum of the 
photon z, respectively; in the absence of rescattering, we have simply = —pi^, see Fig. [2] 
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below. The functions Fe and Fm are given in terms of the proton electric and magnetic form 
factors, via 


Fm{Q^.) = GliQ 


.^2. + QjGim 

4m2 + g2 


(13) 


with 


g|(QD 


GUQl) 

7.78 


1 

(l + Q2/0.71GeV2)^ ’ 


(14) 


in the dipole approximation, where Ge and Gm are the ‘Sachs’ form factors. The modulus 
of the photon virtuality, is given by 


Ql 


qI + 

1-6 


(15) 


i.e. it is cut off at a kinematic minimum gimin = 6 ^^p/(l — 6 )- 

The cross section for the photoproduction of a state V, for the case that the photon is 
emitted from proton i, is then simply given in terms of the photon flux fll 2 p and the 'yp ^ Vp 
subprocess cross section 

(^Ip^pVp = j dA^^(6) (y'^p^vp > (16) 

integrated over the relevant phase space region. As the transverse momentum transferred by 
the photon exchange is typically much smaller than that due to the proton-Pomeron vertex, 
we may safely ignore interference effects, so that the total cross section is simply given by 
summing over i = 1 , 2 , i.e. allowing for the case that the photon is emitted from either 
proton. For two-photon production 77 —)■ X, the corresponding cross section is 


dcr. 


pp^pXp 


dhF. ’ 


dn 


(17) 


77 


where IF™ is the 77 c.m.s. energy. The 77 luminosity is given by 


dL^'T' 


21 F, 


77 


dW^^ dyx 


n{xi) n{x 2 ) , 


(18) 


where yx is the object rapidity and xi ^2 = ^^exp(±|/x), while n{xi) is the photon number 
density: 

n{Xi) ^ . 


(19) 












Figure 2 : Schematic diagrams for the exclusive photoproduction process pp —)■ pVp with (a) 
and without (b) screening corrections included. 


3.2 Soft survival effects 

For photon-mediated processes, survival effects can be included exactly as described in Sec¬ 
tion 12.21 however some additional care is needed. From flTOj) we can see that it is the ampli¬ 
tude for the production process that is the relevant object when including these effects. On 
the other hand, flT^ and flTTl) and the flux flT^ are dehned at the cross section level, with 
the squared amplitude for the photon-initiated subprocesses summed over the (transverse) 
photon polarisations. 

To translate these expressions to the appropriate amplitude level, it is important to in¬ 
clude the photon transverse momentum dependence in the appropriate way, corresponding 
to a correct treatment of the photon polarisation, see [H] . To demonstrate this, we will only 
consider the Fe term in flT^ in what follows, but will comment on the contribution of the 
magnetic form factor at the end. Schematic diagrams for the bare and screened photopro¬ 
duction amplitudes are shown in Fig. |2l with the relevant momenta indicated; for the bare 
amplitude we have qi^ = —pi ^, while for the screened we have = —p'i ^, see Section 12.21 
Using the same decomposition that leads to (jl]), the photoproduction amplitude correspond¬ 
ing to the hgure behaves as 

r(9u) ~ . ( 20 ) 

where is the jp ^ Vp amplitude for a photon of ± helicity, and if the photon is emitted 
from the other proton we simply interchange 1 -f-)- 2. In the bare case {qi± = —pi±) we 
simply square this, and after performing the azimuthal angular integration, the cross terms 
~ Pi p\ vanish and we have 

\T{pi^)\^ pI^(t^p^Vp , ( 21 ) 

where a^p^vp is the subprocess cross section summed over the incoming photon transverse 
polarisations. This is consistent with the Fe term in flT^ and with flT6|) . and indeed a full 
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treatment, keeping all prefactors and expressing the pp —)■ pVp cross section in terms of (| 2 T]) . 
leads to exactly these results, and is the essence of the equivalent photon approximation. 

When calculating the screened amplitude it is crucial to correctly include this explicit 
transverse momentum dependence as in ( 1201 ) . with qi^ = —p'l^ included inside the integral 
dH]); this vector structure of the amplitude can have a signihcant effect on the expected 
survival factor. More precisely for the photoproduction amplitude we take 

T{qi^,q2±) =T' ) + iq\^{A'^ + A )) , (22) 

where T' contains the transverse momentum dependence of the other proton (i.e T' ~ 

([25])), as well as the ^ dependence (and other factors) in flT^ and 
the 7 p c.m.s. energy W^p dependence of the yp —)■ Vp subprocess. 

For the case of two-photon initiated processes, the amplitude can be decomposed precisely 
as in dlj), i.e. 

T{qi^,q2i) ~ ■ ^.2j_)iT++ + T —) + • • • , (23) 

where the T\^\^ are now the 7 (Ai) 7 (A 2 ) —)■ X helicity amplitudes, and we omit the overall 
factors for simplicity. Using this, it can readily be shown that the bare amplitude squared 
reduces to the correct cross section level expressions given in Section [3U1 while in the screened 
case it is again crucial to include this correct vector form of the amplitude, i.e. with the 
qi^ = —p'i^ included inside the integral (|H|). As the relative contributions of the amplitudes 
TajA 2 affect the dependence in fl 2 ^ . the survival factor may depend sensitively on the 
helicity structure of the 77 —)■ X process. For example, in the case of dilepton production 
77 for which the T±± amplitudes vanish for massless leptons, we hnd much less 

suppression than may be naively expected [H] , see Section 14.31 for further discussion. 

Finally, we must also consider the contribution from the magnetic form factor Fm in 
firSj) . While generally suppressed by for larger values of Mx and/or production in the 
forward region, the corresponding value of ^ may not be so small, and the contribution from 
this term may not be negligible. A careful consideration of the derivation of the equivalent 
photon approximation shows that this contribution is generated by a term ~ given by 
the density matrix of the virtual photon (i.e. the proton spin sum) in the cross section. This 
is not proportional to and does not allow a decomposition, at the amplitude level, as in 
fl 20 |) : the Fe contribution on the other hand is given by the term proportional to g(/g//, as 
expected from fl20|) . Therefore, to evaluate the Fm contribution we simply omit any such 
gf^ dependence when calculating the screened amplitude (I 8 |). For the photoproduction case, 
we then add this squared amplitude incoherently to the Fe term, which is calculated as 
described above. For two-photon production, we keep the explicit vector qi^ dependence as 
in fl 2 ^ for the (dominant) Fe{Qi)Fe{Q 2 ) contribution, while for the other terms no explicit 
vector qi^ dependence is included in the amplitude, and the corresponding contributions are 
again squared and added incoherently. 
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Mx(min) 

99 

qq 

bb 

999 

999 

75 

120 

0.073 

0.12 

6.0 

0.14 

150 

4.0 

1.4 X 10-3 

1.7 X 10-3 

0.78 

0.02 

250 

0.13 

5.2 X 10-3 

5.2 X 10-3 

0.018 

5.0 X 10-^ 


Table 1: Parton-level predictions for exclusive two and three jet production cross sections 
(in pb) at the LHC for different cuts on the minimum central system invariant mass Mx 
at \fs = 13 TeV. The jets are required to have transverse momentum > 20 GeV for 
Mx(min) = 75,150 GeV and pj_ > 40 GeV for Mx(min) = 250 GeV and pseudorapidity 
\r]\ < 2.5. The anti-/ct algorithm with jet radius R = 0.6 is used in the three jet case and the 
qq cross sections correspond to one massless quark flavour. Soft survival effects are included 
using model 4 of [50] . 


4 Physics processes 

In the following sections we consider a selection of representative examples of the physics 
processes that are generated by Superchic 2. 

4.1 Exclusive jet production 

Exclusive jet production [5T][52]. in particular of a 2-jet system, has been of great importance 
in testing the underlying perturbative GEP formalism. As discussed in the introduction, it has 
been observed at the Tevatron [T8l[T9] , and there is much potential to measure this process at 
the LHG, in particular with the protons tagged with forward proton spectrometers associated 
with the ATLAS and GMS central detectors, see [15] . Most events with two scattered protons 
and central jets will correspond to central diffractive (GD) jet production, i.e they will not 
be truly exclusive, but will have additional particle production from the Pomeron remnants. 
Exclusive production may be regarded as a particular case of GD jet production with only the 
jets in the hnal state, and no Pomeron remnants. It proceeds through the mechanism shown 
in Fig. [H via the gg —)■ gg, qq and gg —)■ ggg, gqq subprocesses for 2- and 3-jet production, 
respectively. 

Further details about the contributing helicity amplitudes are given in Appendix It 
is found, in particular, that in the case of the gg —)■ qq process, the = 0 amplitude (136]1 
involves a helicity flip along the quark line, and vanishes as the quark mass niq —)■ 0. From 
this fact, we expect a strong suppression in the GEP cross section for quark dijets, relative to 
the gg case, for which the gg —)■ gg amplitudes fl38]) with = 0 incoming gluons display no 
such suppression. In this way the exclusive mode offers the possibility to study almost purely 
(over 99% for typical event selections) gluonic and, crucially, isolated jets [33| produced by 
the collision of a color-singlet gg state, shedding light on their underlying properties (such 
as multiplicity, particle correlations etc) in a well-dehned and comparatively clean exclusive 
environment. In Table [H some representative predictions for exclusive two and three jet 
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production are shown and this gg/qq hierarchy is clear. The invariant mass distributions are 
shown in Fig. [3l with similar results being evident. 



Mx [GeV] 

Figure 3: Parton-level distributions for 2 and 3-jet CEP with respect to the system invariant 
mass Mx at s/s = 13 TeV, using MMHT14 LO PDFs [53]. The hnal-state partons are 
required to he in the pseudorapidity region —2.5 < rj < 2.5 and have transverse momentum 
Px > 20 GeV (leading to minimum invariant masses, Mx, of 40 and 60 GeV in the two 
and three jet cases, respectively), while the three-jet events are dehned using the anti-fc^ 
algorithm with R = 0.6. Distributions are shown for massless quarks, qq, as well as for bb 
production. Soft survival effects are included using model 4 of [50] . 


In the case of three jet production, that is qqg and ggg jets, this suppression in the qq 
exclusive dijet cross section also leads to some interesting predictions [5ll|55]. In particular, 
we expect the behaviour of the qqg amplitude as the radiated gluon becomes soft to be 
governed by the corresponding Born-level, qq, amplitude. As this vanishes for massless 
quarks and = 0 incoming gluons, it is expected to lead to an enhancement of ‘Mercedes- 
like’ conhgurations for the qqg case, where all three partons carry roughly equal energies and 
are well separated. The corresponding three-jet cross sections are also shown in Table [1] while 
the gg dijet cross sections are of order ~ 100 pb, the three-jet ggg cross section are a factor 
of ~ 10 smaller, and the qqg cross section a further order of magnitude smaller again; this 
is due to the specihc colour and spin-dependence of the contributing gg qqg amplitudes, 
which also leads to some suppression in the inclusive case, as well as the additional dynamical 
suppression discussed above. The corresponding invariant mass distributions are shown in 
Fig. [31 

So-called ‘planar radiation zeros’ were shown in [56] to be present in 5-parton QGD 
amplitudes, that is, a complete vanishing of the Born-level amplitudes, independent of the 
particle polarisations, when their momenta he in a plane and satisfy certain additional con¬ 
ditions on their rapidity differences. These were seen in particular to occur in the gg —)■ ggg 
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and, in certain cases, the gg qqg amplitndes, when the initial-state glnons are in a colonr- 
singlet conhguration. This is precisely the sitnation for exclnsive 3-jet prodnction, and so 
it is interesting to examine whether snch zeros may be observable in the CEP process. In 
Fig. m (left) we show parton-level MC predictions for ggg jet prodnction with respect to the 
azimnthal separation Acpij for glnon pairs satisfying the cnt 0.9 < Aij < 1.1, with 



(24) 


where i,j = 1, 2, 3 label the hnal-state glnons, and = yi — yj- It was shown in [56] that a 
zero occnrs when a glnon pairing satishes Aij = 1 and has zero azimnthal angnlar separation 
Acpij = 0. In Fig. m (left), a clear snppression by several orders of magnitnde for lower Acpij, 
which is driven by this zero, is evident; it was shown in [5B] that snch a strong snppression 
is not seen inclnsively. In Fig. 0] (right) we show predictions with respect to the rapidity 
difference \yij — yx\, with the cut Acpij < 10° applied on the angular separation, and with 
the additional requirement that coshAj^ > 4, which helps to isolate the region where A^j = 1 
may be satished. Again a clear dip is seen in the distribution; as before it can be shown 
explicitly that this dip does not occur inclusively, and it is directly driven by the presence of 
a radiation zero, rather than being, say, an artefact of the cut choices. 

While measurements of such distributions can in principle provide quite a clear demon¬ 
stration of these radiation zeros, an important question is whether the expected signal size 
would be large enough. From Table [T] we may expect ~ 500 signal ggg events from 100 pb“^ 
of low pile-up running which can be realistically anticipated at the LHC at s/s = 13 TeV, for 
jet transverse momenta p± > 25 GeV, and rapidities \t]\ < 2.5. By extending out to \t]\ < 5 we 
may increase the event sample by a factor of ~ 1.5, and in addition it should be emphasised 
that the only necessary requirement of the production mechanism for these radiation zeros 
to be present is that the initial-state gluons be in a colour-singlet state. That is, the effect is 
independent of the particle polarisations, and of any additional kinematic effects specihc to 
the pure CEP process, such as the = O’*""'' selection rule. Thus events where one or both 
protons dissociate, but with large rapidity gaps between the dissociation states and the jet 
system, for which the initial-state gluons are also in a colour-singlet conhguration, will be 
expected to contain such zeros. Allowing for these, we may expect a sample of 1000 ggg 
events, while the qqg contribution is expected from Table H to be ~ an order of magnitude 
smaller, and will not be considered in what follows. 

While such a sample may at hrst glance be sufficient to be sensitive to the type of 
dips shown in Fig. 01 some caution is needed. We recall that we are interested in planar 
conhgurations of the jets, where the zeros may occur: after placing a cut of A0jj < 10°(20°) 
the expected sample is reduced to only 55 (142) events. Moreover, it can be shown that the 
zero condition A^j = 1 only has a solution if cosliAjj > 7, corresponding to |Ajjj > 2.6. Such 
a separation in rapidity, which leads to larger invariant masses of the three-jet system, is 
strongly suppressed, and placing such a cut on the original sample of 1000 events reduces 
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Figure 4: Differential cross sections (in pb) at y/s = 13 TeV for exclusive ggg jet production 
at parton level, with respect to (left) the separation in azimuthal angle for gluon pairs 
ij passing the cut 0.9 < Aij <1.1, where Aij is dehned in (124|) . (right) the rapidity differences 
\yi,j ~yx \ for which coshAjj > 4 and A0jj < 10°, where yx is the rapidity of the 3-jet system 
and the ytj is the rapidity of parton i,j, both of which are binned. All partons are required 
to have p± > 20 GeV and \y\ < 5. 


it to only 36. A less restrictive cut may be placed as in Fig. 0] (right), of coshA^ > 4, but 
after combining this with a reasonable cut on Acpij, less than 1 event remains. Clearly a 
measurement of a dip in such a rapidity distribution, even with a more hne-tuned choice of 
cuts, will be highly challenging during low-luminosity running at the LHC. 



0.7 < Aij < 1.3 

0.9 < Aij < 1.1 

Acjij < 60° 

Acjij > 60° 


Acjij < 60° 

Acjij > 60° 


Exclusive 

0.6 

55 

1 .1% 

0.15 

16 

0.9% 

Inclusive 

10.3 

157 

6 .6% 

3.2 

48 

6.7% 


Table 2: Expected number of ggg jet events for which a gluon pairing ij passes cuts on 
the azimuthal separation Acpij and Aij, dehned in the text. Results for both exclusive and 
inclusive parton-level events, using the MMHT14LO PDFs |53] are shown, with an initial 
sample of 1000 events, before cuts, considered in both cases for illustration. 


On the other hand, if the cut 0.7(0.9) < Aij < 1.3(1.1) is placed on the same event sample, 
then 53(16) expected events remain, and even with such a fairly small number of events some 
discrimination may be possible by considering the azimuthal separation observable A0jj as 
in Fig. 0] (left). In Table |2] the expected number of events for which a gluon pairing ij passes 
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such a cut on Aij and with greater and less than 60° is shown; for exclusive production, 
we anticipate a strong suppression in the < 60° region. For comparison results for 

a colour summed inclusive sample of 1000 events are also shown: a strong suppression in 
the ratio of events passing the A0jj < 60° to A0jj > 60° cuts in the exclusive case is seen 
in comparison to the inclusive. Although these results are at LO and parton-level only, 
and may be washed out somewhat in a more realistic treatment, including in particular 
parton shower/higher-order effects as well as background events, this may nonetheless be a 
promising measurement possibility. By comparing such a ratio with the measured sample of 
events dominantly due to CD jet production (for which there is no colour-singlet requirement 
and therefore no observable radiation zeros), a more robust signal of this suppression may 
be observable. 

Finally, it is again worth emphasising that the only requirement on the production mech¬ 
anism for such zeros to occur is that the initial-state gluons be in a colour-singlet state, with 
the particle polarisations playing no role, and so the pure CEP case is not the only possibil¬ 
ity to observe these zeros. As described above, events with proton dissociation may also be 
considered, but it would also be interesting to examine if jet properties such as the colour 
flow and multiplicity might be used to isolate the contribution from a colour-singlet initial 
state, and so investigate these zeros in an inclusive environment. For example, the so-called 
jet pull angle variable has been shown to be useful in identifying jets which originate from a 
colour-singlet gg initial state [57|[58]. 


4.2 Exclusive vector meson photoproduction 


In this section we will consider the photoproduction of vector mesons, focussing on the 
and T(IS') cases; although '0(2*5) production is also included in the MC, it will not 
be considered here. At the LHC, coherent J/il) photoproduction in ultra peripheral p-Pb 
collisions has been measured by ALICE at ^/s]^ = 5.02 TeV [59], and LHCb have made 
increasingly precise measurements oi J/'ip (and '0(25)) photoproduction in pp collisions [91160] 
at \/i = 7 TeV. We will focus here on production (in pp collisions) in the forward region 
relevant to the LHCb acceptance, but will also show some representative results for central 
production. 

Following the notation of Fig. |2l for the yp Vp subprocess cross section we take the 
power-law ht 




N, 


V 


W. 


IP 


IGeV 




by e 


-byql 


(25) 


For the case of J/ip production we take A0, = 3.97 nb and = 0.64, consistently with the 
HERA ht [6T], which Ends = 3.97 ± 0.05 and 5^ = 0.67 ± 0.03; these precise choices will 
be justihed below. For the T(15) we take the values of Ax = 0.12 pb and (5x = 1-6 from [62] . 
although we note that in this case these are quite poorly constrained by the existing HERA 
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data. The slope by is fitted using a Regge-based parameterisation 


+ . ( 26 ) 

with bo = 4.6 GeV~^ and a' = 0.2 GeV~^, consistently with the HERA measurement |63]. In 
the absence of any precise data in the cases of T(IS') and '0(2*5) production, we assume that 
these values are universal. 



2 <p^ < 4.5 

—1 < < 1 


Vi = 7 TeV 

Vi = 13 TeV 

Vi = 13 TeV 

^bare 

359 

511 

333 


278 

406 

291 


0.77 

0.79 

0.87 


Table 3: Gross section predictions (in pb) for exclusive J/'ip —?• photoproduction in 

pp collisions, for different values of the c.m.s. energy ^/s and different cuts on the muon 
pseudorapidities. Results are shown for the ‘bare’ and ‘screened’ cross sections, i.e. excluding 
and including soft survival effects, respectively, and the resulting average suppression due to 
these is also given. 


In Table [3] we show cross sections predictions for J/0 —>■ production at ^/s = 7 and 

13 TeV, with the hnal-state muons restricted to he within the LHGb acceptance {2 < < 

4.5), as well as for central production (—1 < < 1) at ^/s = 13 TeV. The muons are decayed 

including spin correlations, assuming s-channel helicity conservation in the J/'0 production 
subprocess and with the corresponding branching ratio taken from [M]- Predictions are 
shown for demonstration both with and without soft survival effects included, with in the 
latter case model 4 of [50] taken, although the results are in fact almost insensitive to this 
choice. This is to be expected: the main model dependence in the evaluation of the soft 
survival factor lies in the region of small impact parameter b^ Rp, where Rp is the proton 
radius, whereas the peripheral photoproduction process is relatively insensitive to this lower 
bt region. 

The (screened) 7 TeV prediction is in excellent agreement with the LHGb measurement 
of [9] 

(^2 < < 4.5) = 291 ± 7 ± 19 pb , (27) 

where the hrst error is statistical and the second is systematic. However, it is important 
to emphasise that the predicted value depends sensitively on the precise form of the £t in 
(125]) to the 7 P —>■ J/'ipP subprocess cross section, in particular the value of the power 5.^. As 
described above we have chosen a value for this which is at the lower end of the uncertainty 
band of the HERA £t. Taking a larger value will lead to an increase in the predicted cross 
section, with for example 5^ = 0.70, on the upper end of the uncertainty band, giving a 
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dcr/dy^ [nb], ^/s = 7TeV 





Figure 5: Distributions with respect the J/'tp rapidity at i/i = 7 TeV, compared to the 
LHCb data points from [9]. Theory curves corresponding to the ‘bare’ and ‘screened’ cross 
sections, i.e. excluding and including soft survival effects, respectively, are shown, and the 
integrated cross sections are normalised to the data for display purposes. The correlated 
systematic errors are not shown. 


~ 40% larger result, although a more precise evaluation of the uncertainty must account for 
the error on and the anti-correlation between this and 6.^. We therefore choose this value 
to give a good £t to the LHCb data. However, this should be considered as a lower bound 
on the predicted cross sections, due to the low choice of 6^. It is therefore clear from Table [3] 
that without the inclusion of soft survival effects, the LHCb data are in strong tension with 
such a £t to HERA data. 

To examine the influence of survival effects further, we can also consider the distribution 
with respect to the J/ip rapidity, showr0 in Fig. |5l As discussed in Section |2^ the survival 
factor is not constant, and will therefore have an effect on the predicted distributions of 
the hnal-state particles. This is seen clearly in the hgure, with the inclusion of screening 
corrections leading to a steeper fall-off with increasing rapidity. This is to be expected: 
as increases, so does the fractional momentum ^ (for the dominant case 

that the photon is emitted from the proton moving in the positive ^ direction), leading to 
a larger minimum photon see fflSj) . The reaction therefore becomes less peripheral, and 
the survival factor will decrease. This effect is also seen in Table [3l when comparing the 
average survival factor between the central and forward predictions. On the other hand, 
adjusting the input value of 5^ in fl25ll within the range consistent with the HERA data 

^The data points corrected from the fiducial measurement are shown so as to remove the influence of 
the muon cuts, giving a clearer demonstration of the underlying theory; as the correction factors are in fact 
derived in using a previous version of SuperChic, these do not imply any significant model dependence. 
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Figure 6: Distributions with respect the J/iIj transverse momentum at y/s = 7 TeV, cor¬ 
responding to the ‘bare’ and ‘screened’ cross sections for J/il) —)■ production, i.e. 

excluding and including soft survival effects, respectively. The muons are required to have 
pseudorapidity 2 < < 4.5, and the integrated bare cross section is normalised to the 

screened value for display purposes. 


leads to much smaller changes in the predicted distribution. Although the agreement is 
still far from perfedH, the overall trend of the data clearly prefer the screened prediction. 
While this conclusion is only strictly true in the context of the simple power-law HERA £t 
fl25l) . nonetheless this illustrates the importance of a full inclusion of soft survival effects in 
theoretical models such as e.g. [621165] . 

We may also consider the distribution with respect to the J/ip transverse momentum 
This is an important variable in the LHCb measurements [9l|60], for which the selected 
events contain a non-negligible fraction with proton dissociation and/or additional particle 
production that falls outside the LHCb rapidity coverage. To subtract this background the 
measured distribution is fitted by a sum of two exponentials ~ exp(—corresponding 
to the elastic and proton dissociative contributions. The data are well £t by such a parametric 
form, and in [9] LHCb hnd 

6^1 = 5.70 ±0.11 GeV^ (28) 

for the pure elastic CEP contribution, while for the proton dissociative contribution the value 
of the corresponding slope is signihcantly smaller, reflecting the larger average p± in this case. 
Recalling that in CEP the vector sum of the proton transverse momenta is transferred directly 
to the produced object, this fitted value reflects a non-trivial interplay between the elastic 
electromagnetic and Pomeron form factors, given in and fl2^ respectively. While photon 
exchange generally prefers smaller values of the proton p± and so will have a smaller impact 

^However, we note that the correlated systematic errors are not shown in Fig.jSJ 
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on 6ei fl28|) . this contribution cannot necessarily be neglected completely, in particular in the 
forward region where the slope of the Pomeron form factor in fl26H can be quite high, and the 
average photon virtuality is larger. Again, we will expect screening corrections to have some 
influence on this value: in particular, as the expected suppression is larger at higher proton 
p_L, we will expect these to increase 6ei compared to the bare case, see [5] for additional 
discussion. Such an effect, although fairly small, is clearly seen in Fig. [6|, where the J/ijj 
transverse momentum distributions in the screened and bare cases are shown. Performing a 
least-squares £t for < 0.4 GeV^ we find the distributions can be well fitted by a simple 
exponential with slopes 


= 5.0 GeV^ = 5.5 GeV^ (29) 

with a ~ iO.lGeV”^ error due to the uncertainty on the HERA £t [63] in (l26il . and a 
smaller error ~ ±0.02GeV^ due to the fitting procedure. We can see that the bare result 
is inconsistent with the quite precise LHGb measurement (|28|) . but that the introduction 
of survival effects greatly reduces this tension. While the predicted rapidity distributions 
in Fig. [5] and the preference for screening corrections found in that case depend on the 
validity of the power-law fit fl25l) outside the original W^p region of the HERA fit, the 
parameterisation (126]) is grounded in more fundamental principles of Regge theory: the value 
of the slope is driven by the structure of the Pomeron-proton vertex and the slope a' of the 
exchanged Pomeron, while the contribution from the heavy vector boson vertex will be very 
small. This behaviour is therefore expected to be present in more sophisticated models such 
as [62][65]. and thus this result provides a more certain, and less model-dependent indication 
of the importance of a correct, fully differential, inclusion of survival effects. Interestingly, 
it appears that the predicted value may be somewhat lower than the measurement; further 
theoretical investigation of the model dependence of the result, as well as experimentally a 
more precise measurement of &ei, in particular as a function of the J/'ip rapidity, would help 
to clarify this. 



2 < < 4.5 


Vs = 7 + 8 TeV 

y/s = 13 TeV 

yfs = 13 TeV 

T 

^scr. 

0.23 

0.34 

0.29 


Table 4: Gross section predictions (in pb) for exclusive T photoproduction in pp 

collisions, including screening effects, for different values of the c.m.s. energy y/s and different 
cuts on the muon pseudorapidities. The 7 + 8 TeV result is given by the weighted average 
of the predictions corresponding to the relative fractions of integrated luminosity collected 
by the LHGb measurement [66] at these energies. 


Finally, considering the case of exclusive T(IS') —?• production, this has recently 

been measured for the hrst time in hadronic collisions by the LHGb collaboration [66], at ^/s = 
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7 and 8 TeV. As with J/'^ production (|25|) we take a power-law fit to the 7 T cross section, 
but in this case the available HERA data [671168] are much less precise, and the resulting £t 
parameters are not well determined. The default £t in earlier versions of SuperChic, which 
was qualitatively similar to the LO prediction of [65], leads to much too steep an energy 
dependence in comparison to the LHCb data, and a new £t is therefore required. This is 
however complicated by the fact that in pp collisions the photon may be emitted from either 
colliding proton, i.e. we must add the cross sections at energies the 

lower energy W~ contribution, due to the positive power 5 in fl25]l . is expected to be smaller, 
but not necessarily negligible. A simultaneous fit to these contributions must therefore be 
performed, and as the effect of changing 5 typically acts in opposite directions on the 
contributions, this leads to a large amount of cancelation between these, such that a good 
fit to the combined cross section can be achieved for a wide range of 5, in particular for the 
current fairly limited LHCb and HERA data. This highlights a wider issue: when attempting 
to extract information about the jp —)■ Vp cross section from measurements in pp collisions 
some additional assumptions must be made in separating the W~^ and W~ contributions. 
A more robust extraction may be made by increasing the statistics and W range of the 
data, but ultimately a direct and model independent comparison in pp collisions can only 
be performed against lab frame variables, such as the meson rapidity distribution, where 
both and W~ contributions are suitably included in the theory prediction. On the other 
hand, in pA collisions this issue does not arise, as the source of photons can to very good 
approximation be uniquely identified with the heavy ion, due to the enhancement in the 
photon flux. 

Assuming a simple power-law as in fl25]) and fitting to the available data leads to a quite 
low best £t value of h ~ 0.3, but with a sizeable error. Such a low value is disfavoured on 
general grounds [69], from which we would expect a larger 6 than in the case of the lower 
scale J/'0 production process. Moreover, due to the cancellation effects discussed above, it is 
possible to achieve a good £t for much higher values of 6, and instead we set hx = 0.7, with the 
normalisation Nr = 5.7 pb found from the resulting £t, for which we still have x^/d.o.f ~ 1. 
The weighted average of the predicted cross sections within the LHCb acceptance at y/s = 7 
and 8 TeV, corresponding to the relative fractions of integrated luminosity collected by 
LHCb [ 66 | at these energies, is shown in Table H] and is seen to be reasonably consistent 
with the measurement of 0.22 ± 0.07 pb. Predictions for y/s = 13 TeV are also given: these 
contain a ~ 50% uncertainty due to the error in the extracted parameters of the power-law 
fit. 

4.3 Two—photon mediated processes 

In this section we present a very brief selection of results for the two-photon exclusive pro¬ 
duction of lepton (electron and muon) and W boson pair production. In Table |5] we show 
predictions for the muon and W boson pair production cross sections, with and without soft 
survival effects included. In the case of muon pair production we can see that, as expected 
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> 2Mw 

/i'*'/! , <0.1 GeV 

w+w- 

^bare 

6240 

11.2 

3170 

87.5 

<^sc. 

5990 

9.58 

3150 

71.9 

(Si) 

0.96 

0.86 

0.994 

0.82 


Table 5: Cross section predictions (in fb) for exclusive muon and W boson pair production 
at y/s = 13 TeV. The muons are required to have p_L > 5 GeV and \r]\ < 2.5, and are shown 
with and without an additional cut of > 2Mw, while in the W boson case, no cuts 
are imposed. Results are shown for the ‘bare’ and ‘screened’ cross sections, i.e. excluding 
and including soft survival effects, respectively, and the resulting average suppression due to 
these is also given. 


from the discussion in Section [3l the average soft suppression factor is close to unity, due 
to the peripheral two-photon interaction, as well as the vanishing of the T±± amplitudes for 
massless leptons discussed in Section 13.21 However, as seen in the previous section, as the 
system invariant mass increases, we will expect the photon momentum fraction oc Mx 
to increase. This will lead to a higher average photon virtuality, see flT^ . and therefore for 
the average survival factor to be smaller for this less peripheral interaction. We also show 
the prediction for the same muon pair cross section, but subject to the requirement that 
> 2Mw] while the suppression factor is still quite close to unity, it is clearly lower. This 
reduction in the survival factor with Mx is seen more clearly in Fig. [7] where the average 
suppression is shown for muon pair production as a function of the pair invariant mass; a very 
similar result is found in the case of W pair production. We also show in Table [5] the total 
W boson pair production cross section, where the suppression factor is smaller still, due to 
the different helicity structure of the production amplitudes (for which the T±± amplitudes 
are non-vanishing). Finally, the muon pair production cross section, but with the outgoing 
protons required to have transverse momentum < 0.1 GeV is shown: by placing such a 
cut, the reaction is required to be highly peripheral, and it can be seen that the suppression 
factor is extremely close to unity. On the other hand, as discussed in CD], in the case that 
one or both protons dissociate the reaction is generally much less peripheral, and a proper 
inclusion of soft survival effects becomes crucial; this can lead to sizeable deviations in the 
data with respect to the result of e.g. the LPAIR MC [nnn, which does not include these 
effects. 

Recently, the ATLAS collaboration have published a measurement of exclusive and 

production 1731 in normal LHC running conditions, by vetoing on additional charged- 
particle tracks associated with the lepton vertex, and applying further corrections to extract 
the exclusive signal. This is compared to the MC predictions in Table [6l The bare cross 
sections are in both cases too high compared to the data, but a better agreement is achieved 
when survival effects are included. However, interestingly, while there is excellent agreement 
within uncertainties in the electron case, the prediction for the muon cross section lies ~ 3 
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e“'“e 

^bare 

0.795 

0.497 

<^sc. 

0.742 

0.459 

(S.\> 

0.93 

0.92 

ATLAS data [73] 

0.628 ±0.032 ±0.021 

0.428 ±0.035 ±0.018 


Table 6 : Cross section predictions (in pb) for exclusive muon and electron pair production 
at y/s = 7 TeV. The muons (electrons) are required to have p± > 10(12) GeV, and in both 
cases \ri’‘\ < 2.4. Results are shown for the ‘bare’ and ‘screened’ cross sections, i.e. excluding 
and including soft survival effects, respectively, and the resulting average suppression due to 
these is also given. These are compared to the ATLAS data ua. 


a above the data, i.e. a lower value of the average soft suppression appears to be preferred. 
Such a discrepancy may indicate that a further rehnement of the modelling of the opacity 
in the high bt region, to which two-photon induced processes are sensitive, is required, or 
alternatively may be a result of contamination from non-exclusive events due, for example, to 
proton dissociation, although a detailed attempt is made in [73] to subtract this background 
and account for any uncertainty on this in the systematic error on the data. Further measure¬ 
ments, ideally differential in mu, as well as with tagged protons, thus effectively eliminating 
the possibility of proton dissociation, will be of great use in clarifying this issue. It is worth 
emphasising that as the two-photon production process is theoretically so well understood, 
this represents a particularly clean probe of soft survival effects, in particular if the outgoing 
protons are tagged. 

These results highlight the importance of a proper treatment of screening corrections, 
which is still often not included in the literature. In the recent work of m for example, 
where the question of constraining the photon PDF in exclusive and W~^W~ production 
is considered, soft survival effects, which as noted above may be particularly important if 
proton dissociative events are included, are omitteci^. Another important example of this is 
in eg. where an evaluation of the survival factor for two-photon induced processes is given, 
the predictions of which are compared to in the ATLAS data in [73]. While a differential 
treatment of the survival factor is given, and for example the same qualitative decrease with 
Mx as in Fig. [7] is seen, the correct photon q± dependence, described in Section [321 is not 
included in this work; in impact parameter space, only the bt dependence of the photon flux 
flT2|) is included, and not that of the 77 —?• W subprocess. This omits entirely any process- 
dependence in the survival factor, and will not give a reliable estimate for the expected 
suppression. After an explicit calculation, we hnd that including only the bt dependence of 

^Moreover, in the semi-exclusive case it is not the standard photon PDF which enters in the hard cross 
section. Rather, the PDF must be evolved using a modified form of the DGLAP equation in which emission 
in the experimentally relevant rapidity region is forbidden. This will be the subject of a future study m- 
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Figure 7: Average survival factor = dcTscr./clcTbare as a function of the central system 

invariant mass Mx for muon pair production, at ^/s = 14 TeV. The muons are required to 
have p± > 2.5 GeV and \ri\ < 2.5. 


the photon flux, as in izg, tends to underestimate the survival factor by ~ 10 — 20 % for 
lepton pair production within the ATLAS event selection. It should be emphasised that this 
does not correspond to a genuine model dependence: the correct inclusion of the photon q± 
dependence as discussed in Section 13.21 follows simply from the derivation of the equivalent 
photon approximation and cannot be omitted. Indeed, in the case of lepton pair production, 
if we include the correct photon q± dependence at the amplitude level, and instead of using 
the model of [50] we take the simplihed form for the proton opacity used in [76], then the 
predicted survival factor is almost unchanged for these peripheral interactions. 

4.4 Heavy quarkonia production 

The CEP of Xc,b (and rjcfi) quarkonia states has been considered in [5l|6l[77] . These processes 
are implemented in the new MC, and in this section we present updated predictions for a 
short selection of Xc,b cross sections. However, we note that a number of theoretical updates 
and modihcations have been included in comparison to these earlier studies: we summarise 
these below hrst. 

As discussed in these earlier studies, the predicted cross section for Xc CEP has a signif¬ 
icant theoretical uncertainty. One important source of this is the gluon PDF, which at the 
quite low x and values relevant to the process is not well determined; for this reason, 
in 016], consistently with the earlier treatment in [78], the Xc cross section was calculated 
at a lower value of ^/s = 60 GeV, where the gluon is better constrained, and a Regge ex¬ 
trapolation a oc with A = 0.2, was assumed to calculate the cross section at higher 
energies. While the prediction of such an approach was observed to give good agreement 
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with the CDF measurement of XcJ production [79], clearly it depends on the validity of the 
Regge scaling assumption, and in particular on the precise value of A. Such a scaling closely 
corresponds to assuming a typical power-like low-x form xg ~ x~^ in the gluon PDF, which 
appears in the perturbative amplitude ([I|) via (I3|). However, while typically gluon PDFs 
extracted from global hts exhibit to good approximation such a behaviour, the precise value 
of A = d\ogxg/d\ogx, which will depend on the PDF set, may be different from A = 0.2 
and, crucially, will depend on the scale due to DGLAP evolution; such physics will be 
missed by the simple assumption above. 

More physically, we now therefore treat quarkonium production consistently with other 
exclusive processes within the MC, and evaluate the cross section at the appropriate \fs 
value, without any explicit Regge scaling assumption. As the CDF measurement of 77 
CEP [7] strongly favours the results from typical LO gluon PDFs, see e.g. [23], a Ending 
which is further supported by the LHCb measurement of J/V’ pair production at ^/s = 7 and 
8 TeV [22], these should be the appropriate choice here. However, the low x dependence of 
such PDFs leads to a significantly steeper ^/s dependence: for example, with MMHT14LO 
PDFs [53] for a typical given by the average gluon transverse momentum {Q\) ~ 4 GeV^ 
we find A 0.32. The effect of instead taking such a value of A is dramatic, leading to a 
factor of ~ 5 larger cross section at Tevatron energies, and ~ an order of magnitude larger at 
the LHC. The result of explicitly including the corresponding PDF set in ([T]) within the MC 
is consistent with this increase, with the precise energy dependence being somewhat steeper 
still. This leads to a predicted total XcJ cross section at ^/s = 1.96 TeV that is significantly 
larger than the value measured by CDF ra. 

There are however other corrections to the simplified approach of [5l[6] which we may 
consider. In particular, in these works, the Xc,b cross sections are given in terms of the 
99 ^ X vertices 

^ 0 + = (SMj (gi^ g2x) - (gix 92 x ) (gL + 9 IJ- qlj , (30) 

^1+ =-, (31) 

V 2 + = ^^^^’^ (■s(qij/,(q2ja + 2(qi^q2jpi^P2a)e*/'^ , (32) 

where qi± is the transverse momentum of the incoming gluon i, with qfj_ = —qfj_, and 


1 167ra5 / 6 

Crucially, in [5l|6] the vertices fl3nll - fl32l) were calculated in the limit. This 

approximation was justified by the fact that, while a more complete treatment of the particle 
kinematics, as in the /c_L-factorisation approach of the Durham model, includes an important 
part of the higher-order corrections (within collinear factorisation) which generate initial- 
state gluon off-shellness, the qi± ~ region may be particularly sensitive to the exact form 
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of such corrections, for which a full NLO (and beyond) calculation would be necessary to give 
a complete evaluation. Moreover, the derivation of ([T]) is only strictly valid in the 
regime. 

On the other hand, as discussed in EH. see also [HQ], keeping the full qi± dependence in 
99 ^ X amplitudes can have a signihcant effect, and it may still be the case that this gives a 
more reliable estimate of the expected cross sections. Moreover, as we will now show, such a 
treatment is preferred by the existing data on Xc CEP, when the LO PDFs described above 
are used. For this reason, we now choose to include the full gluon qi± dependence everywhere 
in the vertices (IHOjl - ljH^ . The most signihcant effect of this can be seen by observing that 
the term appearing in the denominator of (133|1 

{ 91 ( 12 ) = ^ (Mj - - ql^) , (34) 

which in the on-shell gluon approximation is equal to M^/2, appears to the fourth power in 
the cross section, so that even for relatively small average (q^_|_)/M^, the correction may be 
sizeable. Moreover, in general this ratio is not too small: using the vertex fl30|) in ([ID, for 
example, we have {qjj_)/M^ ~ 0.3. As this will lead to a larger denominator flSTp . the effect 
is to reduce the expected cross section. The precise value of this depends on such factors as 
the object mass, PDF choice, and y/s, but also on the particle spin, through the form of the 
vertex fl30l) - fl3^ . all of which affect (qfj^). Taking y/s = 60 GeV, so as to compare with the 
predictions of [5l|6] , we hnd that including the full qi± dependence in the 99 ^ x amplitudes 
reduces the expected cross sections for central Xcj production by factors of ~ 3, 6 and 2 for 
the J = 0,1 and 2 states, respectively. The reduction in the cross sections is sizeable, and the 
spin dependence is clear, while at larger y/s the higher average gluon leads to a further 
increase in the suppression. For the Xci, the impact is particularly severe: we recall that the 
production amplitude vanishes for on-shell incoming gluons and in the forward proton limit, 
and therefore we would expect the predicted cross section to be particularly sensitive to the 
precise form of the gluon off-shell corrections, as well as the proton spectra. In the case 
of the higher mass Xb the effect is greatly reduced, leading to a factor of ~ 10 - 20 % smaller 
cross section for the Xbo and Xb 2 and a factor of ~ 2 smaller cross section for the Xbi- 

We hnd that this overall reduction in the expected cross sections is largely cancelled by the 
increase due to the use of the LO PDFs at the appropriate y/s value, so that the hnal predicted 
cross sections at the Tevatron are now in good agreement with the CDF measurement ITS] 
of IIct(xcJ —t J/ilix) = 0-97 nb, see Table [7] below. As discussed above, this PDF choice is 
well motivated by the agreement it gives with the CDF measurement of 77 CEP [7| and the 
LHCb measurement oi J/ip pair production [22], as well as being theoretically simpler and 
arguably more justihed then invoking additional Regge scaling arguments. The full inclusion 
of off-shell gluon effects then appears to be required to give good agreement with the CDF 
Xc data ua when such LO PDFs are used. It should however be emphasised that while the 
agreement with these data is equally good between this updated approach and that of [5l|6], 
the predictions for other observables are quite distinct. This different treatment of the gluon 
off-shellness affects the Xci and Xc 2 to Xco ratios, and the predicted energy dependence of 
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the Xc cross sections is steeper, leading to increased rates at the LHC. In addition, as the 
higher mass Xb states are mnch less sensitive to off-shell glnon effects we predict larger cross 
sections for these, see Table m below. 

A farther change compared to the earlier works relates to the normalisation of the deriva¬ 
tive of the wave fnnction at the origin (j)'p{0) in fl3^ . As in |6] this is normalized in terms of 
the total Xco width 

rtot(Xo) ~ r(xo ^ 99 ) = 96-^|0p(O)|^ . (35) 

This implicitly inclndes higher order contribntions to the 99 ^ x vertex, nnder the assnmp- 
tion that these are the same as in the x ^ 99 case. While in [B] an additional snbtraction 
was made so that these corrections were only inclnded for the spin-zero Xo^ by dividing the 
extracted valne of |0p(O)p for the Xi ,2 by a K-factor, taken to be 1.5, we now choose to apply 
this correction nniformly; in this case the ratio of cross sections between the spin states cor¬ 
respond to the pnrely LO results, consistently with the rest of the calculation. This change 
leads to a factor of 1.5 smaller Xco cross section compared to previous results. In the case of 
the Xb, the measurement of |8T] suggests that earlier estimates for the total width, such as 
that taken in di, underestimate the true value, and prefers rtot(xw) ~ 0.8 MeV, excluding 
any NLO K-factor. We take this value here: as in m this was taken as the value including 
a K-factor (again assumed to be 1.5), the Xbo predictions are unchanged, while the Xb(i,2) 
predictions are now a factor of 1.5 larger. 

Finally, while in [5l[6] an additional ‘non-perturbative’ contribution was considered in the 
case of Xc production, this is no longer included. This does not rule out the possibility of 
important non-perturbative effects for the relatively low mass Xc states; rather, we choose to 
consistently only consider the perturbative mechanism in the MC. The influence of possible 
non-perturbative effects can then be seen as deviations from these predictions, for example 
in the ratio of cross sections between the different spin states. In addition, as discussed 
above, we recall that no enhanced survival effects are considered here: these may reduce the 
expected cross sections somewhat, but will leave the ratios Xi/Xo etc almost unchanged. 


cr X Br 

XcO 

Xcl 

Xc2 

S XcJ 

XcO 

y.xcj 

CTEQ6L1 

0.29 

0.12 

0.22 

0.63 

0.46 

MMHT14LO 

0.38 

0.11 

0.29 

0.78 

0.49 

NNPDF3.0 

0.51 

0.15 

0.41 

1.1 

0.48 


Table 7: Differential cross section, da/dy^, predictions (in nb) at y^ = 0 for exclusive Xcj —t 
production at \/i = 1.96 TeV, including screening effects. The sum over the three spin 
states, and the ratio of the Xco to the total cross sections are also given. Results are shown for 
three choices of LO PDF: CTEQ6L1 [82], MMHT14 [53], NNPDF 3.0 (^^(Ml) = 0.130) [83], 
and with model 1 of [50] for the soft survival factor. 
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cr X Br 

XcO 

Xcl 

Xc2 

Xcl 

XcO 

Xc2 

XcO 

CTEQ6L1 

19 

8.3 

15 

0.44 

0.79 

MMHT14LO 

51 

17 

41 

0.33 

0.80 

NNPDF3.0 

23 

7.1 

19 

0.31 

0.83 

LHCb data 

9.3 ±4.5 

16.4 ± 7.1 

28.0 ± 12.3 

1.8 ±0.7 

3.0 ±0.9 


Table 8: Cross section predictions (in pb) for exclusive Xcj production at 

y/s = 7 TeV, including screening effects, with the cross section ratios between the different 
spin states shown for clarity. Results are shown for three choices of LO PDF: CTEQ6L1 |82j . 
MMHT14 [53], NNPDF 3.0 (as(M|) = 0.130) [H3|, with model 1 of [50] for the soft survival 
factor, and with the Xc decay products in the rapidity region 2 < rj < 4.5. The preliminary 
LHCb measurement [8] is shown for comparison, with statistical, systematic and luminosity 
errors added in quadrature; for the cross section ratios only the statistical errors are included. 


Turning now to the MC predictions, in Table [7] we show the predicted differential cross 
sections at = 0 for Xcj production at the Tevatron for a range of PDF sets, and with model 
1 of the soft survival factor from [50] . This can be compared to the previous predictions of [6] . 
As described above, the steeper energy dependence due to the LO gluon PDF sets and the 
suppression due to the new treatment of off-shell gluon effects are found to largely cancel, 
such that the total XcJ cross section is similar in size, and in all cases consistent with the 
CDF measurement na. On the other hand the ratios Xc(i, 2 )/Xco are different. From the 
discussion above, in the case of the Xa/Xco ratio, we expect the updated off-shell gluon 
effects to reduce this by a factor of ~ 3, while the Xc 2 /Xco ratio should be increased by a 
factor of ~ 2, and in both cases the new treatment of the NLO K-factor for the Xco will 
increase these by a factor of 1.5. We therefore expect the Xci/Xco ratio to be decreased by 
a factor of ~ 2 and the Xc 2 /Xco ratio to be increased by a factor of ~ 3: this is indeed 
consistent with the results of Table [71 It should be emphasised that as the changes due to 
the updated treatment of the off-shell gluon effects are theoretically motived, and essential 
to reproduce the CDF Xc data ca with the LO PDFs sets preferred by measurements of 
exclusive processes at the Tevatron and LHC [2^17^ (as well as themselves being a consistent 
choice for the LO theoretical calculation), these may be expected to better describe the data; 
the difference due to the treatment of the NLO K-factor, on the other hand, which enters 
at the level of the uncertainty due to higher-order corrections, simply corresponds to one 
possible prescription, with everything treated at LO, while the approach [B], where some 
approximate K-factors are included, may give a better estimate. Finally, we note that the 
ratio Xco/ 'YhXcJi also shown in Tabled is reasonably consistent with the more recent CDF 
limit for this of < 50%, calculated from the corresponding limits on Xco production in the 
7r'''7r“ and K^K~ channels [84] . 

In Table |8] we show predictions for Xc J/V'T —t CEP at y/s = 7 TeV for 
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cr X Hr 

XcO 

Xcl 

Xc2 

Xcl 

XcO 

Xc2 

XcO 

CTEQ6L1 

30 

14 

26 

0.46 

0.87 

MMHT14LO 

77 

24 

68 

0.31 

0.88 

NNPDF3.0 

39 

11 

35 

0.28 

0.90 


Table 9: Cross section predictions (in pb) for exclusive XcJ production at 

^/s = 13 TeV, including screening effects, with the cross section ratios between the different 
spin states shown for clarity. Results are shown for three choices of LO PDF: CTEQ6L1 |82j . 
MMHT14 [33], NNPDF 3.0 = 0.130) [M], with model 1 of [SU] for the soft survival 

factor, and with the Xc decay products in the rapidity region 2 < rj < 4.5. 


the same event selection as the preliminary LHCb measurement |8] of this process. While 
the default predictions with the MMHT14LO PDFs are somewhat too large, the agreement 
is better, although not perfect, with the CTEQ6L and NNPDF3.0 PDF sets. These are 
on the other hand only the predictions for the central PDFs: for the MMHT14LO set, for 
example, the PDF uncertainty is found to be ± ~ 50%. The Xc{i, 2 )/Xco ratios also exhibit 
some generally small PDF dependence, which is larger in the xa case. These ratios are 
qualitatively consistent with the predictions, and it should be emphasised that this is in itself 
a highly non-trivial result: for example, for the Xc 2 without the effect of the Jz = 0 selection 
rule we would expect the corresponding ratio to be ~ 2 orders of magnitude higher, as is 
observed inclusively. However there remains some indication of tension with the LHCb data 
in both cross section ratios, which, as discussed in detail in |6], may be due to the influence of 
proton dissociation on the selected events, but may also indicate that further theoretical work 
is needed to better model the CEP of these relatively low mass states. Further measurements 
of this, with the new HERSCHEL detectors [I6] at LHCb allowing a much greater rejection 
of non-exclusive events, will greatly clarify this question; the corresponding cross section 
predictions at y/s = 13 TeV are shown in Table HI 

A further observable, which is not expected to suffer from the same uncertainties as in 
the Xc case, is the CEP of the heavier Xb states: for example, as discussed above, the effect of 
adjusting the treatment of gluon off-shellness is much smaller, due to the larger meson mass, 
and more generally the production process is safely in the perturbative regime. In Table [10] 
we show predictions for Xbj ^ Ty —)■ production at ^/s = 13 TeV, within the LHCb 

acceptance. Due to the larger mass, the predicted cross sections are signihcantly smaller, 
although still experimentally realistic. The expected Xbi/Xbo ratio is much smaller than in the 
Xc case, due to the higher meson mass, while the Xb 2 /Xbo ratio is somewhat smaller, mainly 
due to the differing branching ratios for the particle decays. We note that the predicted cross 
sections are signihcantly larger than the estimates of [6]: the SuperChic vl.47 MC predicts 
for example a Xbo cross section of 8 fb within the same acceptance. This sizeable difference 
is primarily due to the steeper energy dependence induced by the LO PDF sets. 

Finally, it is interesting to observe the dependence of the average survival factor on the 
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cr X Br 

XbO 

Xbl 

Xb2 

X 61 

XbO 

Xb2 

XbO 

CTEQ6L1 

73 

5.3 

27 

0.073 

0.37 

MMHT14LO 

110 

6.6 

43 

0.060 

0.39 

NNPDF3.0 

62 

3.5 

24 

0.056 

0.39 


Table 10: Cross section predictions (in fb) for exclusive XbJ —)■ T 7 —)■ production 

at y/s = 13 TeV, including screening effects, with the cross section ratios between the dif¬ 
ferent spin states shown. Results are shown for three choices of LO PDF: CTEQ6L1 |82j . 
MMHT14 [33], NNPDF 3.0 {as{M'^) = 0.130) [M], with model 1 of [SU] for the soft survival 
factor, and with the xt decay products in the rapidity region 2 < 7 < 4.5. 


spin on the produced states. This is shown in Table fTTl for Xc production at y/s = 1.96 and 
7 TeV, and for two different soft models in [30]: the average suppression is much weaker for 
the higher spin states, due to the more peripheral nature of the interaction vertices fl3T]) and 
fl32|) . while the stronger overall reduction as the collider energy is increased, and the variation 
between the model choices, is clear. 



Model 1 

Model 4 

XcO 

Xcl 

Xc2 

XcO 

Xcl 

Xc2 

(>S'eik), Tevatron 

0.027 

0.067 

0.060 

0.033 

0.088 

0.077 

(Sk), LHCb 

0.014 

0.029 

0.036 

0.019 

0.041 

0.050 


Table 11: Average survival factors for XcJ production corresponding to Tables [7] (‘Tevatron’) 
and [ 8 ] (‘LHCb’), for two different model choices of the soft survival factor dehned in |30] . 
Values are shown for MMHT14 LO [53] PDFs, but the results for other choices are similar. 


4.5 Generated processes: summary 

In the coming years the LHC will take increasingly precise data at unprecedented energies. 
This presents the possibility for a wide programme of exclusive measurements, building on 
those already performed during Run-I [12]. For this reason a wide range of processes, some 
of which have been discussed in detail above, are included in the MC. We summarise these, 
and give some motivations for further experimental measurements, below (see [311H11231IH3] for 
further details of the processes generated): 

• Exclusive jet and vector meson production, discussed further in Sections 14.11 and 14.21 
respectively. 
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• XcJ production, discussed in Section WM by both the J/V’7 ^ind two-body 

decay channels, including full spin correlations. As discussed in Section 031 the existing 
preliminary LHCb measurement [ 8 ] of exclusive Xc J/'07 production, while quali¬ 
tatively agreeing with the Durham model expectations, shows some possible tension 
with the theory, in particular with an apparent excess of Xc 2 events visible. Further 
measurements of this process, with the new HERSCHEL detectors [16] at LHCb allow¬ 
ing a much greater rejection of non-exclusive events, will greatly clarify this question. 
Observables such as two body decays, where the Xc 2 is no longer enhanced by the J/ijJX 
branching (and the Xci will be absent entirely), as well as more differential observables 
such as the Xc(i, 2 )/Xco ratio as a function of p±(x) will allow a closer comparison to 
theory. 

• Light meson pair (vrTr, KK, , (p4>) production. Measurements of the meson p± 

distribution are of particular importance, as these are sensitive to the transition from 
the ‘non-perturbative’ low p± region where a Regge-based model may be more appro¬ 
priate to the higher p± region where the perturbative approach described in [87] should 
become relevant. Where such a transition occurs is currently a very uncertain ques¬ 
tion. However, we note that the perturbative approach makes very distinct predictions 
about, for example, the cross section for flavour-non-singlet meson pair production, 
compared to Regge-based expectations. 

• 77 production. As well as being of interest in its own right, this relatively clean 
process can serve as a useful benchmark to compare with other CEP measurements: 
by considering cross section ratios with respect to 77 CEP, uncertainties due to the 
gluon PDFs and survival factor largely cancel. In addition, as discussed in [ 88 ] the 
light-by-light scattering process has yet to be seen experimentally, and may represent 
an observable at the LHC, or in a future circular collider; the photon-induced 77 —)■ 77 
process is therefore also included. 

• Double J/^ and t/’(2S') quarkonia production. In [22] LHCb have reported the ob¬ 
servation of double J/'ip and J/ijj'ip{2S) production, in broad agreement with theory 
expectations in [23]. More precise measurements, in particular in the case of 'ip{2S)J/'ip 
and double 'ip{2S) production (so far not seen), as well as more differential measure¬ 
ments of the double J/'ip cross sections, for example with respect to the meson rapidity 
separation and would provide a useful test of the perturbative approach. As dis¬ 
cussed above, a measurement of the ratio to other exclusive processes such as 77 would 
also greatly reduce the potentially large theoretical uncertainties. In addition, although 
not currently included in the MC, the cross section for double pc production could be of 
the same size or even larger than for the J/'0, due to the particular Feynman diagrams 
which enter in this case. 

^Such a model is not included in SuperChic: see [5^[5S] for MC implementations. 
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• Xb production, by both the Ty —)■ /i+p “7 and two-body decay channels, including full 
spin correlations. The higher mass compared to the Xc case means that the process 
should be safely perturbative, so that we should expect good agreement between the 
measured Xb 2 /Xbo ratio and theory, while the Xbi should be practically absent. 

• Pseudoscalar rjc^ meson production. These are expected to be strongly suppressed by 
the = O’*" selection rule. A measurement of these processes would therefore provide 
a direct test of this; for the r/c, the predicted cross sections are large enough that this 
could represent a realistic observable, although it may be challenging to identify an 
experimentally viable decay channel. 

• W^W~ and production (Z = e, p,r), via two-photon collisions, see Section 031 
This purely QED cross section is known theoretically to very high precision, and in the 
case of lepton pair production, if cuts are place to restrict the proton (or equivalently 

to low values the survival factor may be essentially omitted, so that such a 
process could be used as a luminosity calibration [IH]. Alternatively, without placing 
such cuts, this well understood process can serve as a particularly clean probe of soft 
survival effects, in particular in the presence of tagged protons. Exclusive W^W~ 
production is of particular interest as a probe of potential anomalous gauge couplings, 
see [89] - l9T| . All photon-induced processes are also available for initial states. 

• SM Higgs boson production. The CEP of a SM Higgs, in precisely the mass region 
observed by ATLAS and CMS [92],[93], has received a lot of attention. As discussed 
further in [3] , the exclusive mode is very well suited to probe crucial identihcation issues 
such as the bb coupling and the CP-parity of this object. 

• The photoproduction of p(770) and 0(1020) mesons, via the tt+tt" and K~^K~ decay 
channels, respectively. As for other photoproduction processes, the simple hts fl25l) and 
fl26p for the W dependence of the cross section and slope parameter b are taken. The 
fit of [SI] to the effective Pomeron trajectories is used, and the normalisations are set 
using data from [95l|96]. These give 6 = 0.19 (0.16) and a' = 0.125 (0.158) GeV“^, 
with a = 11.4 (0.96) pb at W^p = 75 (70) GeV and b = 11.1 (7.3) GeV“^ at W^p = 84 
(70) GeV, for the p (0). 

In addition, although not currently included in the MG, the measurement of ‘exotic’ 
quarkonia in the exclusive channel may provide further information about these states. For 
example, the X(3915), for which there has been recent discussion as to whether this can be 
interpreted as a Xco(2P) state [971[98], could be searched for exclusively in the DD channel. 
The Z(3930), currently associated with the Xc 2 i‘^P), rnay also be searched for via the DD 
mode; as mentioned in m, the expected cross section for such an excited state should be 
similar to the Xc 2 {^P)- In addition, as discussed further in |3], a clear observation of exclu¬ 
sive X(3872) production can provide important constraints on its nature. Other possibilities 
include the J/00 narrow resonances, such as the X(4350), clearly seen by Belle [99], for 
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which O’*"’*' or 2^^^ assignments, as well as tetraqnark or more conventional charmonia inter¬ 
pretations, remain possible. Snch a resonance peak could be searched for in the exclusive 
channel, as part of more general resonance searches in for instance the ip'cf) channel. A further 
process not currently included in the MC is exclusive DD production, for which quite large 
cross sections are expected, and is sensitive to the cc —)■ DD transition with no additional 
particle production. This latter process is the subject of a future work [75], including a full 
MC implementation. 

5 Super Chic 2: code and availability 

SuperChic 2 is a Fortran based Monte Carlo that can generate the processes described in 
Section IT31 with and without soft survival effects. User-defined distributions may be output, 
as well as unweighted events in the HEPEVT and Les Houches formats. The code and a user 
manual can be found at http://projects.hepforge.org/superchic. 


6 Summary and outlook 

In this paper we have presented results from a new Monte Carlo for central exclusive pro¬ 
duction (CEP): SuperChic 2. In a CEP process, an object X is produced, separated by 
two large rapidity gaps from intact outgoing protons, with no additional hadronic activity. 
Theoretically, the study of CEP requires the development of a framework which is quite 
different from that used to describe the inclusive processes more commonly considered at 
hadron colliders, while experimentally it results in a very clean signal (in the absence of pile 
up), and the outgoing protons can be reconstructed by proton tagging detectors, situated far 
down the beam line. 

SuperChic 2 builds on the earlier SuperChic event generator, but includes a range of 
theoretical updates and improvements. Most significant of these relates to the soft survival 
factor, representing the probability of no additional underlying event activity, which would 
spoil the exclusivity of the final state. In previous MC implementations, this has simply 
been treated as a constant quantity. However, this is not in general true, and we have 
demonstrated in this paper a number of cases where such an assumption breaks down. For 
the first time, SuperChic 2 includes a fully differential treatment of the soft survival factor, 
and we have seen that this can lead to some distinct, and model-dependent, predictions for 
the corresponding particle distributions. 

In addition, a much wider range of processes has been included, with exclusive 2 and 3 
jet production being a particularly relevant addition; this is the first inclusion of the latter 
process in a publicly available MC. Other production processes currently included are: vector 
meson (p, (/>, J/V', V'(2*S'), T(IS')) photoproduction, Xc,b and rjcfi quarkonia, photon pairs, light 
meson pairs, heavy '0(25')) quarkonia pairs, SM Higgs boson and two-photon initiated 

lepton, W boson and photon pair production (with both pp and initial states). 
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We have described the theoretical approach implemented in the SuperChic 2 MC, and 
have considered in more detail a selection of the processes listed above: in Section ed] we have 
considered 2 and 3-jet production, and discussed the possibility for observing the ‘radiation 
zeros’ described in [56] in the 3 jet case; in Section 021 we have considered J/'0 and T(IS') 
photoproduction, showing how a correct inclusion of soft survival effects leads to important 
differences in the predicted meson rapidity and transverse momentum distributions; in Sec¬ 
tion 14.31 we have considered two-photon initiated muon and W boson pair production, and 
shown how the predicted survival factor depends on the system mass Mx] in Section 03] we 
have described an updated approach to modelling quarkonium CEP, and presented results 
for Xc,b production. 

In light of the measurement possibilities for exclusive processes during Run-II of the 
LHC [12], the study of CEP processes is very topical: exclusive events may be measured 
with both protons tagged using the approved and installed AFP [13] and CT-PPS [14] 
forward proton spectrometers, associated with the ATLAS and CMS central detectors, re¬ 
spectively, see also [TB], as well as using rapidity gap vetoes to select dominantly exclusive 
event samples, for example at LHCb, where the relatively low instantaneous luminosity and 
wide rapidity coverage allowed by the newly installed HERSCHEL forward detectors [16] are 
highly favourable. It is our aim that the new SuperChic 2 MC will play an important role 
in this promising physics programme. 
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A Jet production: helicity amplitudes 


In this appendix we give some further details about the contributing amplitudes in the 
exclusive jet cross section. The helicity amplitudes relevant to quark dijet production are 


■^±±,hh 

■^±T,hh 


6‘ 


cd 


lOvra, 




2AT, (1 - /12 cos2 e) Mx 
J'”’* Svra, 


{I3h ± , 
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while for the gluon amplitudes we have 
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-^±=F,+- 
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- 1 (1 - cos2 9) ’ 




CD 


Nc 

-r- StOOs 

m-1 


1 ± cos 9 \ 

1 =F cos 0 y ’ 


(38) 

(39) 


for gluons of ‘±’ helicity and quarks of helicity h/2, while c, d are the outgoing quark color 
labels, and (3 = {1 — Arn^/. As discussed in [23], these amplitudes are shown for 
the azimuthal angle 0 = 0 for the outgoing particle momenta. In general, some helicity 
amplitudes have an overall 0-dependent phase, which while having no effect in the spin- 
summed inclusive cross section, must be included when using the decomposition (jl|) to give 
the correct result in the exclusive case. 

The amplitudes relevant to 3-jet production are all ‘maximally helicity violating’ (MHV) [lOUj . 
and can therefore be written down quite simply, see [SB] for further details. 
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